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ABSTRACT 

We present CO(l-O) maps of 12 warm H 2 -selected Hickson Compact Groups (HCGs), covering 14 individually 
imaged warm H 2 bright galaxies, with the Combined Array for Research in Millimeter Astronomy. We found a 
variety of molecular gas distributions within the HCGs, including regularly rotating disks, bars, rings, tidal tails, 
and possibly nuclear outflows, though the molecular gas morphologies are more consistent with spirals and 
early-type galaxies than mergers and interacting systems. Our CO-imaged HCG galaxies, when plotted on the 
Kennicutt-Schmidt relation, shows star formation suppression of (^) = 10±5, distributed bimodally, with flve 
objects exhibiting suppressions of ^ > 10 and depletion timescales > 10 Gyr. This star formation inefficiency 
is also seen in the efficiency per freefall time of|Krumholz et al.|(|2012]i. We investigate the gas-to-dust ratios of 
these galaxies to determine if an incorrect Lco-A/(H 2 ) conversion caused the apparent suppression and And 
that HCGs have normal gas-to-dust ratios. It is likely that the cause of the apparent suppression in these objects 
is associated with shocks injecting turbulence into the molecular gas, supported by the fact that the required 
turbulent injection luminosity is consistent with the bright H 2 luminosity reported by |Cluver et al. ( |2013| ). 
Galaxies with high star formation suppression ^10) also appear to be those in the most advanced stages 
of transition across both optical and infrared color space. This supports the idea that at least some galaxies in 
HCGs are transitioning objects, where a disruption of the existing molecular gas in the system suppresses star 
formation by inhibiting the molecular gas from collapsing and forming stars efficiently. These observations, 
combined with recent work on poststarburst galaxies with molecular reservoirs, indicates that galaxies do not 
need to expel their molecular reservoirs prior to quenching star formation and transitioning from blue spirals 
to red early-type galaxies. This may imply that star formation quenching can occur without the need to starve 
a galaxy of cold gas first. 

Keywords: galaxies: evolution — galaxies: star formation — galaxies: kinematics and dynamics 


1. INTRODUCTION 


The present-day galaxy population has a bimodal distribu¬ 
tion, comprised of blue spiral galaxies and red elliptical and 
lenticular galaxies ( |Tinsley|1978[|Strateva et al.|2001[|Baldry| 
et al.|2004 ), with a dear th of ^axies at intermediate optical 


colors ( [Bell et al.|[2003| ). Their rarit y suggests that galaxies 
transition rapidly in colors. However, Schawinski et al.|pQ14| ) 
showed that selecting transitioning objects based on color 
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leads to an overestimation of morphologically transforming 
objects, as the majority of green valley objects are not those 
undergoing the rapid morphological transformation from spi¬ 
ral to elliptical, but were mainly spiral galax ies that were un¬ 
dergoi ng secular evolution. More recently, [Smethurst et al. 
( |2015| ) supported this picture of spiral galaxies transitioning 
at intermediate rates in the green valley, but also showed that 
early type galaxies transition more rapidly. 

Many early transitioning scenarios posited that the majority 
of molecular gas in galaxies is depleted prior to the quench- 
ing of star formation ( [Sanders & Mirabel] 1996t[Hopkins et aT 
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Table 1 

HCG CO(l-O) galaxy properties 


HCG 

Name 

Principal 

Name 

RA 

(J2000) 

Dec 

(J2000) 

Dist.^ 

(Mpc) 

Morph. 

t-type 

^1.4GHz 

(mJy) 

[Cii]? 

MOHEG^ 

Optical 

Sequence 

Spitzer 

gap 

WISE 

IRTZ 

Optical'^ 

AGN 

co^ 

Morph. 

HCG 25b 

PGC012539 

03 2045.41 

-010240.9 

85.8 

0.9±1.3 

5.6“ 

/ 

/ 

red 

/ 

/ 

/ 

D 

HCG 40c 

PGC027508 

09 38 53.61 

-04 5136.6 

94.1 

2.3±2.1 

10.3^ 



green 



/ 

B+R 

HCG 47a 

UGC05644 

1025 46.26 

+13 43 00.7 

141 

2.4±1.2 

13.1“ 



blue 



/ 

R,S 

HCG 55c 

PGC035573 

113205.69 

+7048 38.7 

222 

1.3±1.3 

- 

/ 


green 




D 

HCG 57a 

NGC3753 

11 37 53.90 

+2158 53.0 

127t 

2.1±0.6 

3.8“ 

/ 

/ 

red 

/ 

/ 

/ 

M 

HCG57d 

NGC3754 

1137 54.92 

+2159 07.8 

127t 

3.4±1.1 

- 

/ 


blue 




R 

HCG 68a 

NGC5353 

13 53 26.69 

+4016 58.9 

34.6 

-2.0±0.7 

40.5“ 


/ 

green 



/ 

D 

HCG 68b 

NGC5354 

13 53 26.70 

+4018 09.9 

38.1 

-2.1±0.7 

8.0“ 

/ 

/ 

- 

- 

- 

- 

- 

HCG 79a 

NGC6027A 

15 5911.14 

+204517.5 

64.5 

0.3±2.1 

9.3“ 

/ 


red 


/ 

/ 

D 

HCG 82b 

NGC6163 

16 28 27.91 

+325047.0 

148 

-1.7±1.1 

- 


/ 

red 

/ 

/ 

/ 

M 

HCG 91a 

NGC7214 

22 0907.68 

-27 48 34.1 

92.6 

4.5±0.7 

29.2“ 

/ 


blue 



/ 

S 

HCG 95c 

PGC071077 

23 19 31.09 

+09 3010.7 

158 

9.0±2.0 

4.10^ 

/ 

/ 

green 

/ 


/ 

M 

HCG 96a 

NGC7674 

23 27 56.72 

+084644.5 

116t 

3.8±0.6 

221.0“ 

/ 


blue 



/ 

B+R,S 

HCG 96c 

PGC071505 

23 27 58.78 

+084658.1 

116t 

5.5±4.7 

0.85^ 

/ 


green 




D 

HCG 100a 

NGC7803 

0001 19.97 

+13 0640.5 

69.5 

O.lil.O 

12.3“ 

/ 


blue 



/ 

D 


^ Luminosity distance determined using the Nearby Extragalactic Database (NED) 

^Distance determined by the more massive of the gr oup members 

Morphological t-type deter mined by HyperLEDA ^Makarov et al.|2014l _ _ _ _ _ 

Fia GHz continuum from (a) |Condon et al.|(|1998), (b)|Menon & Hickson|ri985j, (c)|Brown et ^(2011) , and (d) |Becker et al.|{T9^ . 
o MOHEG definition based on H2/ 7.7jUm > 0.04 (Ogle et al.|2UU6[|Cluver et al.|2013) i, all except 57d are detected m H 2 . 

^ Optical AGN classification from |Martmez et al. \2010) , counting both transition objects (TO), LINERs, and Seyferts 

^Morphological cla ss of the molecular gas: D = disk, R = ring, M = mildly disrupted, B+R = bar/ring, and S = spiral, based on the classification scheme from 
lAlatalo et alljM^ 


|2006| ), through both supernov a and active galactic nucleus 
(AGN) feedback mechanisms ( [Springel et al.|[2Q05) ). How¬ 
ever, recent observations have started to question whether the 
two-stage scenario consisting of (1) the expulsion of star¬ 
forming gas followed by (2) the cessation of star formation 
(SF) is the exclusive evolutionary picture. Studies now show 


cone et al.|2012[[MT4^|Alatalo|2015| ), and tidal disruption and 
harassment through group inte ractions ([Hickson et al.|jl992 


Zabludoff & Mulchaey||1998^ [Walker et al.||2010[ [Bitsakis 
et al.|2010 |2014j ). Given their low galaxy velocity d ispersion 


and h igh aensity, as well as their short-lived nature ( [Hickson 
1982[ ), compact groups serve as an ideal environment in which 


of molecular gas ([French et al.[[2015 

[Rowlands et aL[[2015) 

and are dustier than normal galaxies ( 

Yesuf et aL|2014| 

). This 


shows that removing all star-forming material before transi- 
tioning is not a requirement in a galaxy’s transformation from 
blue to red. 

With the advent of the Wi de-field Infrared Survey Explorer 
(WISE; [Wright et al.[[201Q} mission, evidence has accumu¬ 


lated that mid-infrared (IR) colors can also be used to iden- 


tify phases of transitioning galaxies ^Ko et al.f|2013j) . [Alatalo 
[et al.[([2014b[) showed the existence of a prominent bifurcation 


to study galaxy transformation. 

Hickson compact groups (HCGs) are defined as “small, rel¬ 
atively isolated systems of typic ally four or five gala xies in 
close proximity to one another” ( [Hickson|1982[ 1997| ). They 
tend to have a high fraction of early-type galaxies (E/SO), evi¬ 
dence of tidal interac tions, and high density structure with low 
velocity dispersion dHickson 1982[ 1997[) and deficiencies 
in H i compared with isolated ga laxies ( [Verdes-Montenegro[ 
et al.|2Q0f|[Borthakur et al.|2010[ ). Compact groups appear to 


[et al.[ ( [MT4^ showed the existence of a prominent bifurcation 
between star-forming spiral galaxies and quiescent early-type 
galaxies in the WISE [4.6]-[12]/im colors, deeming this the 
“Infrared Transition Zone” (IRTZ). Objects within the IRTZ 
have red optical colors (also described in [Ko et al.[2Q13| ), sug¬ 
gesting that galaxies traverse the optical green valley before 
the IRTZ. Assuming that [4.6]-[12]/im WISE colors trace the 
interstellar medium (ISM) within the ga laxy and optical col - 
ors trace the mean stellar population age ( [Donoso et al.[2012] ), 
which further supports the idea that a non-negligible number 
of galaxies quench SF before shedding their ISMs. 

Many plausible mechanisms have been introd uced to ex¬ 
plain this transition, including ma jor mergers ( [Toomre &| 
Toomrep^972[ Springel et al. 2005); experiencing ram pres 


go through an evolutionary phase that can be tra ced by neutral 
gas depletion ( [Verdes-Montenegro et al.[[2001[ ). In the later 
stages of depletion, H l is found less i n the galaxies and mor e 
in the intragalactic medium (IGM); ( Borthakur et al.[[2010[ ), 
with a rise in the fraction of groups containin g extended 
group-wide X-ray emission ( Ponman et al.|199^ . However, 
the origin of the exten ded X-ray emission is sT iII unclear for 
HCGs in general (see Rasmussen et al.[[2008[ ), with strong 
starburst winds being the p artial cause in at least one system 
([O’Sullivan et al. |2014a|b[ ). The fraction of galaxy types also 


evolves in concert with the neutral gas depletion, with spiral- 
rich groups a t early times and elliptical-rich groups later in 
the sequence ( [Bitsakis et al.|2010][2011[[2QT4| ). 

It also appears that the galaxies within the HCG transform 


sure Str ipping, and strangulation when falling into a ga laxy 
cluster ( [Bek£ et 'aL|[2002t [Blanton & Moustaka^[2009| and 


rapidly. No Spitze r color bimodality was reported by [Lacy 


2009[ 2013[); minor mergers ( 

Qu et al. 2010[ Eliche-Moral 

etal.[2012 

[Alatalo et al.[2014c 

mm 

); AGN feedback (Hop- 

kins et al. 


Fischer et a 


IFeruglio et al. 

POTOT 

Sturm et a 


[Alatalo et al.|2011 

[Aalto et al.[[20] 



et al. (2004 2007]) between early-type and late-type galaxies 
m a sample of field galaxies, but [Johnson et al.[ ( 2007[ ) doc¬ 


umented a marked Spitzer IR “gap” m compact group galax 
ies, with very few galaxies obs erved between the st ar-forming 
cloud and the quiescent cloud. [Walker et al.[ ( [2010] ) suggested 
that compact group galaxies within this “gap” were l ikely in 
the midst of a rapid morphological transformation. [Cluver 
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Figure 1. A comparison between the total CO luminosity measured by 
CARMA and IRAM. Errors for CARMA include both the RMS and the 20% 
millimeter flux calibration uncertainty. For IRAM, the range of values con¬ 
sists of the luminosity measured in the central pointing as a lower limit and 
the extrapolated total luminosity as the higher estimate. Overall, the CO lu¬ 
minosities derived from CARMA and IRAM are in good agreement. 


et al. ( |2013| ) showed that the “gap” galaxies in compact groups 
tend to have warm hydrog en emission traced by the Spitzer In¬ 
frared Spectrograph (IRS, |Houck et al.|2004| ) that is enhanced 
beyond the level that photon-dominated regions powered by 
star formation alone could explain, t ermed Molecular Hydro- 
gen Emission Galaxies (MOHEGs; |Ogle et aL]|2007| ). |Clu-| 
|ver et al.| ( |2013| ) suggested that this enhanced H 2 emission 
might be energized by shocks caused by collisions with the 
clumpy intragroup medium. MOHEG “gap” HCGs therefore 
represent ideal laboratories to test the interplay between rapid 
galaxy transition, SF quenching, and the disruption of the star¬ 
forming fuel and ISM. 

Evidence that star formation might be taking place ineffi¬ 
ciently in transitioning galaxies has started to mount, includ 


ing in radio galaxies ( [Nesvadba et al.|[20T0 


2015|), AGN-driv en molecular outfiow host NGC 1266 ( |Alat 


Guillard et al. 


alo et al.llMBbl ), as well as in HCGs. Both the shock in 
Stephan’s Quintet ([Appleton et al.|20Q6}|Guillard et al.|2009| 


2012a[ |Konstantopoulos et al.||26l4|) and the HCG galaxy 
HCG57a ( [Alatalo et al.||2014aj i seem to exhibit suppressed 
star formation. In these cases, the authors suggested that the 
injection of turbulence might be causing this star formation 
inhibition, either from collisions with the radio jets or AGN 
outfiows, or the collisions and interactions within the group 
environment. 

We present new Combined Array from Research in Mil¬ 
limeter Astronomy (CARMA|^CO(1-0) observations of 12 
HCGs, including multiple galaxies considered “gap” galax¬ 
ies and MOHEGs (shown in Table [T]). In ^ we describe the 
sample selection and observations from CARMA, including 
reduction and analysis methods. In we present the molec¬ 
ular properties of the sample, incl uding their posit ion on the 
Kennicutt-Schmidt (K-S) relation ( |Kennicutt|1998| ). In Q we 
discuss these results in the context of transitioning galaxies. 


^ http://www.mmarray.org 


In ^ we summarize our results. We comment and present 
the maps of individual galaxies in the Appendix. The cosmo- 
logical parameters H p = 70 km s“^, = 0.3 and = 0.7 

( Spergel et al.|2007 ) are used throughout. 


2. OBSERVATIONS AND ANALYSIS 
2.1. Sample Selection 

While single dish observations are able to provide informa¬ 
tion about the CO luminosity, and possibly some information 
about the molecular gas kinematics (such as line-width and 
line-shape), it is unable to provide spatial information about 
the molecular gas. Interferometry is able to provide informa¬ 
tion about the extent and distribution of the molecular gas, 
allowing for direct comparisons between the CO and stellar 
mass, or star formation. For this reason, we chose to use 
CARMA to follow up a subset o f the CQ-bright , warm H 2 - 
bright HCG galaxies presented in |Lisenfeld et al.| ( [2014| ). 

The galaxies chosen for the new CARMA imaging were 
drawn from the HCGs that were detected in warm H 2 using 
the Spi tzer IRS[^ (usually with boosted THa,wa rm/f^ 7 . 7 umPAH 
ratios; Ogle et ar||2007[ [Cluver et al |2013 ) as well as de¬ 
tected m CO(l-O) with t helnstitut de Radioastr onomie Mil- 
limetrique (IRAM) 30m ( [Lisenfeld et al.||2Q14|). On ly 10% 
of HCG galaxies are MOHEGs ( [Cluver et al.||2013] ). Erom 
this subset of CO-bright, H 2 -strong galaxies, we chose a sub¬ 
set that maximized the overlap with Herschel [C ll]-detected 
objects (Appleton et al., in preparation). These objects also 
tended to show early-typ e morphologies based on their t- 
types from HyperLED^^ (Makarov et al. 20141. 11 HCGs 
altogether were chosen for the CARMA sample, including 
8 HCGs with corresponding [C ll] and [O l] Herschel obser¬ 
vations. HCG 96, although not observed through the Spitzer 
HCG program, had data in both the CARMA and Herschel 
archives, and was therefore added to our sample. All sample 
galaxies and their properties are listed in Table 

2.2. CARMA 

The HCGs were observed with CARMA, an interferomet¬ 
ric array of 15 radio dishes (6 x 10.4 m and 9x6Tm) lo cated 
in the Eastern Sierras in California ( |Bock et ST||2006| ). Ob¬ 
servations were taken over the course of three semesters be¬ 
tween 12 Mar 2013 - 16 Jun 2014. One, HCG 96, was taken 
from the archive. Thus, we observed a total of 12 HCGs (and 
15 individual galaxies, including HCG 68b, detected only in 
3mm continuum^ The observing strategy and data reduction 
were perfo rmed in a manner id entical to that of the ATLAS 
galaxies in | Alatalo et al. { 201^ . Tablej^presents the semester, 
gain calibrator, bandpass calibrator, total hours on source, 
and beam full-width at half-maximum (EWHM) for each of 
the HCGs observed with CARMA in the D-array (with base¬ 
lines between 11-150m). Galaxies with new CARMA ob¬ 
servations had a correlator configuration of 8 x 500 MHz win¬ 
dow in each sidebandjj with the CO(l-O) line utilizing the 
highest resolution 500 MHz mode (with channel resolution of 
^15 km s“^). This meant there was a sufficiently large band¬ 
width to measure 3mm continuum, which was successfully 


^ One galaxy, HCG 57d, was not included in the original Spitzer IRS foot¬ 
print. 

^ http://leda.univ-lyonl.fr/ 

^ HCG 40e was also tentatively detected, but below a signal-to-noise ratio 
of 3, and only a small subset of the velocity structure was recovered. 

^ The observations of HCG 96 were taken with a correlator conflguration 
of 5 X 250 and 3 x 500 MHz per sideband, totalling 5500 MHz of bandwidth. 
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Table 2 

CARMA CO(l-O) Observing Parameters 


Name 

Semester 

Total 

Hours'’ 

Gain 

calibrator 

^maj ^ ^min 
(") 

Kelvin 
per Jy 

AV^ 
(km/ s) 

HCG 25 

2014a 

12.14 

0339-017 

6.1 X 3.6 

4.2734 

21 

HCG 40 

2014a 

17.67 

0825+031 

5.9 X 3.9 

4.0316 

21 

HCG 47 

2014a 

12.41 

0956+252 

4.9 X 3.4 

5.4678 

32 

HCG 55 

2013b 

7.64 

1048+717 

7.4 X 7.1 

1.7767 

33 

HCG 57 

2013a 

8.61 

1224+213 

4.6 X 3.3 

6.0303 

31 

HCG 68 

2013b 

16.24 

1310+323 

5.1 X 3.7 

4.9006 

30 

HCG 79 

2013b 

9.75 

1613+342 

7.4 X 6.4 

1.9520 

31 

HCG 82 

2014a 

10.83 

1635+381 

4.3 X 3.6 

5.9389 

21 

HCG 91 

2013b 

8.06 

2258-279 

8.5 X 4.1 

2.6671 

21 

HCG 95 

2013a 

13.26 

3C454.3 

3.9 X 3.0 

7.9045 

21 

HCG96t 

2010b 

7.33 

-0010+109 

2.5 X 2.4 

15.0203 

10 

HCG 100 

2013b 

9.75 

3C454.3 

4.6 X 3.6 

5.5428 

21 


^Channel width 

'’Total time on-source time 

t Archival CARMA data 


Table 3 

CARMA Derived Properties 


Name 

r’sys 
(km/ s) 

Zco 

Vel. Range 
(km/ s) 

rpt 

^100 GHz 

(mJy) 

RMS^ 
(mJy / bm) 

Flux* 

(Jy km/ s) 

Leo 

(lO^Le) 

(1O®M0) 

Area 

(□") (kpc2) 

(Mo/pc2) 

25b 

6371 

0.02125 

6037-6684 

< 0.525 

4.83 

21.16 ±2.37 

1.83 ±0.20 

1.60 ±0.18 

191 

30.4 

52.7 ±5.9 

40c 

6419 

0.02141 

6085-6753 

< 0.634 

6.69 

93.32 ±2.97 

9.70 ±0.31 

8.49 ± 0.27 

563 

107.7 

78.9 ±2.5 

47a 

9637 

0.03215 

9429-9844 

< 0.307 

3.25 

64.56 ± 1.89 

14.91 ±0.44 

13.06 ±0.38 

992 

408.4 

32.0 ±0.9 

55c 

15760 

0.05257 

15394-16125 

< 0.595 

3.97 

10.09 ± 1.33 

5.67 ±0.74 

4.96 ± 0.65 

247 

233.1 

21.3 ±2.8 

57a 

8723 

0.02910 

8214-9231 

< 0.391 

3.85 

59.16 ±3.94 

11.12 ±0.74 

9.74 ± 0.65 

487 

164.6 

59.1 ±3.9 

57d 

8944 

0.02983 

8753-9134 

< 0.391 

3.85 

7.50 ±0.33 

1.41 ±0.06 

1.23 ±0.05 

176 

59.3 

20.8 ± 0.9 

68a 

2249 

0.00750 

1457-3040 

6.687 ±0.198 

5.07 

23.65 ±2.16 

0.34 ±0.03 

0.30 ±0.03 

152 

4.2 

71.1 ±6.5 

68bt 

- 

- 

- 

1.597 ±0.243 

- 

- 

- 

- 

- 

- 

- 

79a 

4149 

0.01384 

3841-4457 

0.756 ±0.172 

5.02 

25.42 ± 1.91 

1.25 ±0.09 

1.10 ±0.08 

382 

35.4 

31.0 ±2.3 

82b 

10415 

0.03474 

10181-10652 

0.483 ±0.126 

3.69 

12.92 ± 1.28 

3.28 ±0.32 

2.87 ±0.28 

99 

44.5 

64.6 ± 6.4 

91a 

6857 

0.02287 

6543-7171 

< 0.991 

9.53 

173.54 ± 4.80 

17.44 ±0.48 

15.28 ±0.42 

2050 

377.4 

40.5 ± 1.1 

95c 

11540 

0.03849 

11270-11831 

0.710 ±0.169 

4.15 

24.74 ± 1.66 

7.13 ±0.48 

6.25 ± 0.42 

173 

87.3 

71.6 ±4.8 

96a 

8638 

0.02881 

8482-8793 

3.911 ±0.383 

7.74 

121.75 ± 1.59 

19.09 ±0.25 

16.72 ±0.22 

533 

150.5 

111.1 ± 1.5 

96c 

8809 

0.02938 

8592-9025 

< 1.150 

7.74 

9.61 ±0.13 

1.59 ±0.02 

1.39 ±0.02 

18 

5.3 

260.3 ±3.5 

100a 

5220 

0.01741 

4888-5551 

0.429 ± 0.090 

4.32 

36.99 ± 2.34 

2.11 ±0.13 

1.84 ±0.12 

343 

36.3 

50.7 ±3.2 


^Does not include 20% absolute flux calibration uncertainty 
^RMS noise per channel _ 

^Does not include 30% conversion uncertainty jBolatto et al.|2013} i n ago, assuming ATdq =2x 10^° cm“^ (K km s“^)“^ 

^Detected only in 3mm continuum; in the CARMA observations, see |Lisenfeld et al. 1(2014) for its detected CO properties based on IRAM data 


detected in 7/15 sources (listed in Table [^. Continuum con- 
tribution s were subtracted in uv-spsicc u sing the MIRIAD task 
uvlin ( |Sault et al.|1995| ), as detailed in |Alatalo et al.| ( |MT3] ). 

The resulting channel maps an d moment maps wer e also 
constructed in identical fashion to I Alatalo et al.|(|2013b. Fig¬ 
ures [CTHCTJ] showcase the CARMA data for each individual 
HCG, including the channel maps, integrated intensity (mo- 
mentO), and mean velocity (moment 1) maps, as well as in¬ 
tegrated spectra and position-velocity diagrams (PVD). The 
PVD was constructed by creating a slice in the velocity cube 
at a certain position angle (shown as a dashed line on the mo- 
mentO figure), and integrating across a slice in space. 

The integrated spectrum was constructed by using the mo- 
mentO map to create a clip-mask and integrating the flux 
within the momentO-defined (unmasked) aperture. This was 
done separately for each galaxy. The root mean square (RMS) 
noise was then taken by calculating the standard deviation 
of all pixels in the cube that were outside of the momentO- 
aperture per channel and is listed in Table An additional 
noise correction of 30% was also added in quadrature to the 


RMS noise to accou nt for the oversampling of the maps (see: 
Alatalo et al.]2015a for details). The RMS noise per channel 


for the integrated spectrum was then calculated by multiply¬ 
ing the RMS of the entire data cube by the square root of the 
total number of beams in the momentO-aperture. 

To calculate an integrated line fiux for each galaxy, we de¬ 
termined the extent of emission within the cube (shaded blue 
in the integrated spectrum) and summed across them. The 
line fiux RMS was then calculated by multiplying the RMS 
per channel of the spectrum by the velocity width (listed in 
Table and the square root of the total number of channels 
integrated to derive the line fiux. 

Table[3]lists the properties derived from the imaging data of 
the individually detected galaxies, including the RMS noise in 
the channel maps, the detected continuum levels, and the total 
detected CO line fiuxes. The spatial extent of the molecu¬ 
lar gas was determined by summing the total number of un¬ 
masked pixels (the momentO-aperture) in the moment maps. 
CO luminosities (shown in Table were calculated using the 
luminosity distance to the source (listed in Table [^. The H 2 
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Figure 2. (Left): SFR and gas s urface density in our HCG sample are shown in comp arison on the Ke nnicutt-Schmidt relation jKennicutt| 199 8) . T he SFRs were 


taken froml Bitsakis et al.|(|2014), normal ized to a Salpeter Initial Mass Funct ion (IMF;|Salpeter|1955|. The H CGs are compared t o the Milky Way (|Yuse f-Zad^ 
et al. 12009 , normal galaxies or LIRGs (|Kennicutt|1998||Fisher et ^ CO-imaged early-type galaxies JDavis et al.|2014), high redshift objects dGenzel 

et al.|2UlU|, and radio galaxies Jogle et al.|2UlU||Alatalo k al.|2Ulbbf, all renormalized to a Salpeter IMF The solid black line represents the Kennicutt-ScTmudt 
relation, and dashed lines represent levels ot 1/10 and 1/lUU suppress' ion (below) and enhancement (above). The shapes and colors of the HCG points are based 
on whether they meet the MOHEG criterion from Ogle et aL| d2007b, with red squares for HCGs with L^.r^mPAK > 0-04, considered confirmed 

ihold. Ihe S 


MOHEGs, and blue circles representing those below this th reshold. Ihe SER error bars include the 10% erro rs quoted in|Bit sakis et a r]j2014) , and the 0.38 dex 
scatter in MAGPHYS-derived SERs from|Lanz et al.|(|2013), added in quadrature, (^ ght): The SE laws from|Kruniholz et aLW'ZOl'Z} are plotted for nearby and 
high- 2 : galaxies (black dots and gray dia monds, respectiv ely; |Krumholz et al.|2012|and references therein), and early-typ e galaxies (indigo square s; [Davis et al.| 

' eter|1955|. The axes ot panel b were chosen to match the axes shown inj&umholz et al.U2012F Dotted lines 

HCGs found to be suppressed in panel a are also suppressed m panel b. 


represent suppressions and enhancements ot 10 and 100 
mass was determined using the LC 0 -MQI 2 ) relation: 


M(H2) = 1.05 X 10' 


^ ScqAvDI 
1 “h f^sys/ 


Mo, 


( 1 ) 


which assumes Xqo = 2 x 10^^ cm ^ (K km s ^ (the 


mean conversion factor presented in [Bolatto et al.||201^ . 
Sco^'^ is the CO(l-O) flux (in Jy km s“^), is the lu- 
minosity distance (in Mpc), Vsys is the optically-defined sys¬ 
temic velocity (in km s“^), and c is the speed of light (in 
km s“^). The time variability of fiux calibrators at 3mm con¬ 
tinuum adds an additional « 20% uncertainty to the CO(l-O) 
fiux, and the Lqq-M(112) conversion also carries a 30% un¬ 


certainty ( [Bolatto et al.||2013] ), imposing an additional 35% 
absolute flux uncertainty to the measured CO masse^ 

2.3. Comparison of the CARMA and IRAM data 

Thirteen of the 14 CQ-image d HCG galaxies presen ted in 
this paper were also observed by jLisenfdd et al.| ( |2014| ) using 
the IRAM 30m (only HCG 96c was not observed). Figure H 
compares the CO(l-O) luminosities derived from both sets m 
data. The measured CO(l-O) luminosities from CARMA and 
the IRAM 30m are in good agreement, confirming that the 
CARMA observations do not resolve out substantial fiux, and 
are therefore a reasonable representation of the molecular gas 
in these systems. HCG 57d, for which CARMA detected only 

^ Errors reported in Table f^do not include the 35% uncertainty, though it 
is included in error bars in alHgures. 


half of the fiux detected with the IRAM 30m (p ossibly due to 
resolving out fiux, and sensitivity), is an outlier (jAlatalo et al. 
[20T4|). 


3. RESULTS 


3.1. Molecular gas morphologies of the HCG galaxies 
Following the morphological classification used by [Alat- 


alo et al. (2013 ) for ATLAS^^ early-type galaxies, we have 
morphologically classified the molecular gas in all of the CO- 
imaged HCGs (listed in Table detailed discussion can be 
found in Appendix |W. We classify each galaxy as being ei¬ 
ther a disk (D), spiral (S), bar-fring (B-fR), ring (R), mildly 
disrupted (M), or a combination of multiple distinct classifi¬ 
cations. Details of the morp hological classifications can be 
found in jAlatalo et al. ( 2013| ). The morphologies seen in our 
HCGs tends to be a niix of regular rotation and dynamically 
excited structures, with a lack of strongly disrupted objects, 
which are quite prevalent in Ultralum inous Infrared Galaxies 
(ULIRGs; [Sanders & Mirabel| 1 9961) and interacting galaxies 
( Wilson et al.|2008[ ). Overall though, our CO-imaged galax- 


ies have app earances comparab le to field galax ies, including 
both spirals ([Heifer et al.|2003[) and early-types (jAlatalo et al.[ 

[ 20 T 3 ] ). 

3.2. The Kennicutt-Schmidt relation in warm, H 2 -bright 
HCG galaxies 

Figure [^displays the molecular gas surface density of each 
HCG compared with its star formation rate (SFR) surface den¬ 
sity (calculated assuming the SF and gas are co-spatial), and 
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Figure 3. The molecular gas-to-dust ratios of the HCGs in this sample, us- 
ing the molecular gas mass determined in this work, a nd the dust masses from 
jBitsakis et al.|(|2014|, with an uncertainty of 0.5 dex dda Cunha et al.|2010J. 
1 he dust mass of HCG 55c is an upper limit. Red and blue dotted lines repre- 
sent factors of two of the average gas-to-dust ratio. The colors of the points 
represent the suppression of SF in the molecular gas in each galaxy. The gas- 
to-dust ratios of the HCGs appear consistent across a large range in molecular 
gas masses, as well as values of .5^. 


using the SFR derived in Bitsakis et al.| (|2Q14| ). SFRs cal¬ 
culated by MAGPH YS (|da Cunha et al.|2008| ) seem to have a 

, wh 


scatter of 0.38 dex ( |Lanz et al.|2013| ), which we have included 
in our SF r ate uncertainty. To t est whether the SF histories 
assumed in [Bitsakis et al. ( |2014| ) could significantly alter our 
conclusions, we input varying SF histories (including continu- 
ous and truncated m odels) and derived SF rates using CIGALE 
( jCiesla et~aL]|2015| ). The SF rates derived from CIGAL E did 


not vary by more than the assumed scatter reported in Lanz 
|etaLl ( [20T3] ). 

^S FR fof the HCGs is comp ared with User of the Milky 
Way ( Yusef-Zadeh et al. |2009 ), normal galaxi es ([Kennicutt 
Fisher et al.|2013| ), high redshift gala xies (|Genzel et"ar 



ratios from Cluver et al. ((2013), with MOHEGs represented 
by red squares, and non-MOHEGs represented by blue cir¬ 
cles. Overall, most HCG galax ies fall within the scatter of the 
K-S relation ( |Kennicutt|1998| ), but HCG galaxies as a group 
are more likely to be found in the lower half of the scatter than 
the upper half (although a larger sample will be required to 
determine whether this is statistically signific ant). This result 
is consistent with the single-dish results from (Lisenf eld et ^ 


( |2014 ), and there are some objects (in particular HCG 25b, 
40c, 57a, 79a, and 82b) that exhibit a much more substantial 
scatter off of the K-S relation. We define the degree of SE 
suppression, as the ratio between the expected SER sur¬ 
face density (Esfr) from the measu red molecular ga s surface 
density Emoi using the K-S relation (|Kennicutt|1998|), and the 
observed Esfr- Values of ^ for each HCG are listed in Ta- 
blelll 

The molecular gas depletion (the molecular gas mass di¬ 
vided by the star formation rate) timescales associated with 
the HCG galaxies are also shown in Table|^ As expected, ob¬ 


jects with normal ^ appear to have depletion timescales sim- 
ilar to those found i n normal, star-forming g alaxies (^ 1 Gyr; 
Bigiel et al.||2008| [Saintonge et al.||^llb| ), and early-type 
galaxies ( (Davis et al.||2014) . The objects with ^ >10 tend 
to have depletion timescales that are of order a Hubble time 
or longer, with HCG 57a and HCG 82b having depletion 
timescales of > 30 Gyr. These depletion timescales suggest 
that these systems will nearly indefinitely contain molecu¬ 
lar gas, if the suppressions present are continuous (rather 
than discrete and sporadic) processes. Given the timescale 
for galaxies within Hickson Compact groups to transform 
(- IGyr; (Walker et al.][20TQ| (Bitsakis et aLl[MTQ| [MTT] ), it 
is likely that the suppression is also a transient phenomenon. 

In order to compare to early-type galaxies, another sup¬ 
pressed set of objects, we have calculated ^ for CQ(l-O) 
imaged ATLAS^^ galaxies ( Young et al.|2008 [Crocker et aL 


(2011 [(Alatalo et al!](2013[ (Davis et al.((26l4( ). We calculated 

the total SERs of these objects using the (uncorrected) WISE 
22/im emissionj^ assuming that the majority of the 22/im 
emission originat es from SE in most C O-bright objects with¬ 
out strong AGNs ( (Calzetti et al.|2QQ7] Nyland et al., submit- 
ted). In the CO-im aged subsample ATLAS^^ galaxies from 
avis et al.((2Q14(), we find ^^ 2.6, consistent with the pre- 
ctions of (Martig et~ar] ( [201 3( ) that the bulges in early-type 
galaxies can stabilize molecular disks, creating suppressions 
of 2-5. 

Eor the HCGs in this work, we find the mean 
(^hcg) ~ 10±5. There also appears to be a marked bi¬ 
modality seen between galaxies forming stars at normal effi¬ 
ciency, and a few that are extremely inefficient (with ^ > 10), 
including HCG 25b, 40c, 57a, 79a, and 82b. Given the com¬ 
pact natur e of the mo lecular gas and strong X-ray emission in 
HCG 68a ( Liu|20lf ), it is quite possible that a non-negligible 
fraction or its rar-infrared emission originates from the AGN, 
rather than from SE, similar to NGC 1266 ( [Alatalo et al.( 
2015b| ). The contributions of AGNs to the spectral energy 


distributions of individual galaxies is discussed in Appendix 

s 

We have also plotted the SE relation following the proce¬ 
dure in [Krumholz et al.[ ( [2012[ ). These authors claim that the 
SE rate in normal (non-LIRG) molecular configurations is de¬ 
pendent on the local freefall time of each individual molecular 
cloud, conve rting about 1% of the m olecular mass to stars per 
freefall time. (Krumholz et al.(([2012[) argue that since the ma¬ 
jority of stars form in giant molecular clouds (GMCs), their 
individual properties are the dominant determinant of the star 
formation efficiency in the molecular gas of a galaxy, and thus 
the local GMC freefall time i s more importan t than the global 
dynamical time (as is used by (Kennicutt 1998[ ). Assuming that 
the G MC distribution within sw athes of galaxies is fairly con¬ 
sistent, (Krumholz et al.( ( [2012( ) use this framework to univer¬ 
salize star formatio n laws, across norm al, star-forming galax¬ 
ies and interactions. (Davis et ar] ( (2Q14| ) have also been able to 
show that this also applies successfully to early-type galaxies, 
while (UFomo et al^ ( (2015( ) have shown that the GMCs within 
example early-type galaxy, NGC 4526, matc hes the distribu¬ 
tions o f late-type local galaxies. We use the Krumholz et al. 
( (2012| ) relation to see if our star formation suppression can be 
reconciled by examining the GMC-scale efficiencies, rather 
than the global scales. This framework was able to link Milky 

^ The scatter associated with the ATg-band based fact or used to correct the 
22[im emission for aged stars from|Davis et al.U2014J is almost as large as 
the observed inefficiency. 
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Table 4 

Star formation suppression values 


Name 

-j-Cl 

(Myr) 

c^h 
^ KS 

^ eff 

+d 

^dep 

(Gyr) 

^norm SF 

(km/s) 

^ci=l 
(km/s) 

(10^5 ergs) 

T ^ 

^turb 

(10^^ erg/s) 


25b 

10.3 

12.2 

9.1 

10.0 

91 

51 

7.2 

31 

0.34 

40c 

11.4 

18.2 

10.4 

12.7 

104 

62 

57 

303 

0.36 

47a 

9.1 

3.8 

3.9 

3.8 

39 

40 

35 

120 

0.50 

55c 

8.2 

3.2 

4.3 

3.8 

43 

32 

9.0 

25 

2.15 

57a 

10.6 

38.9 

26.8 

30.4 

268 

54 

49 

226 

0.82 

57d 

8.2 

0.6 

0.8 

0.7 

8 

32 

2.2 

6.0 

- 

68a 

11.1 

1.0 

0.6 

0.8 

6 

59 

1.8 

9.0 

0.75 

79a 

9.0 

15.4 

16.2 

15.7 

162 

39 

2.9 

9.6 

1.32 

82b 

10.8 

42.2 

27.4 

31.9 

274 

56 

16 

76 

0.34 

91a 

9.7 

1.4 

1.3 

1.3 

13 

44 

52 

200 

0.33 

95c 

11.1 

3.6 

2.2 

2.6 

22 

59 

37 

193 

0.57 

96a 

12.4 

1.6 

0.7 

1.0 

7 

74 

158 

997 

0.14 

96c 

15.4 

3.4 

0.9 

1.5 

9 

113 

31 

297 

0.22 

100a 

10.2 

1.1 

0.8 

0.9 

8 

50 

7.9 

34 

0.95 


Derived SF suppression parameters for the molecular gas in HCGs including (a) th e free-free timescale usi ng Equation2, assuming that (jgas = 10km s ^,(b) 
5^ (the deviation from the K-S relation, (c) the star formation suppression from the|Krumholz et aLU2012| law using fff (and assuming (jgas = 10 km s“^) and 
(d) the depletion timescale (Mmoi/SFR) for each of the CO-imaged HCG galaxy, (e) llie neces sary molecular velocity dispersion (crgas,norm sf) of the giant 
molecular clouds within the galaxies in order for each galaxy to have normal SE efficiency from|Krumholz et al.U2012y (/) The required global molecular gas 
velocity dispersion to stabilize the molecular gas against collapse (Toomre Q > 1). (^f) The corresponding total turbulent energy injection required to stabilize 
the molecul ar disk, {h) The turb ulent luminosity (the total turbulent energy injected per rotational period), (i) The ratio of the total warm H 2 luminosity from 
Spitzer IRS ( jCluver et al.|2013) to the required turbulent luminosity. 


Way, local group galaxies, starburst, and high-redshift galaxy 
scales. The free- fall time for each HCG was calculated using 
Equation 4 froni |Krumholz et^ ( 2012| ): 


= 


rl/4 


78 


( 2 ) 


where G is the gravitational constant, Hgmc is the average 
surface density of the giant molec ular clouds within th e sys¬ 
tem (estimated to be 85 Mq pc“^; Bqlatto et al. [2008 1 . Ugas 
for each object is presented in Table]^ and a is the velocity 
dispersion of the molecular gas within the individual GMCs, 
which we have assu med to be 10 km s~^ (c onsistent with the 


Krumholz et al. 


8 km s ^ assumed in 
pear to vary by more than a factor ot two m disk galax ies ( |Dib 


2012|). (Jgas does not ap 


|et al.|[2Q0^ Walter et al. j2008[ [Chung et al.||2009} , though 
(Jgas ^10 km s“^ might be an undere stimate given the dis - 
rupted kinematics seen in some HCGs ( jAlatalo et al.|2Q14a| . 
Figure shows that some HCGs do not appear to fall on 
the typical efficiency (molecular mass converted into stars per 
free-fall time), e^O.Ol relation with the other galaxies, and 
those that are suppressed in the original K-S plot (Fig[^) con¬ 
tinue to appear suppressed in SF efficiency space as well. The 
least efficient galaxy, HC G 82b, is 27 times less ef ficient than 
the mean efficiency from [Krumholz et al.j ( [2012] ). A larger 
sample of objects with suppressions at this level is needed to 
be able to constrain a duty cycle and determine how long ob¬ 
jects are observed with these large depletion times. 


3.3. Gas-to-dust ratios in H 2 -bright HCG galaxies 

Could the observed suppression in these HCGs b e the re- 
sult of using an incorrect Lco-M(H 2 ) conversion? [Downes 
[& Solomon[ ( |1998[ ) were able to show that ULIRGs did not 
follow the standard Milky Way relation. The molecular gas 
in ULIRGs is not distributed in discrete GMCs, instead be¬ 
ing a more continuous distribution of molecular gas, and the 
velocity widths associated with the additional gas motions 
boosted the CO luminosity per unit gas mass. Anomalous 
Lco-M(Y{ 2 ) conversion factors have also been identified in 


other MOHEGs, including NGC4258 ( [Ogle et al.|2Q14[ ) and 
3C 293 ( [Lanz et al.|2015[ ), and thus is an important factor to 
check in our sample. 

In order to determine whether our observations overesti¬ 
mate the molecular gas mass in the suppressed HCG galax¬ 
ies, we turn to the gas-to-dust ratio as a test. Figure shows 
the gas-to-dust ratios of th e CO-imaged HCG sa mp^ using 
dust masses calculated by [Bitsakis et al.[ ( [2014[) using full 
UV-to -sub-mm SED fitting m MAGPHYS ( da Cunha et al.[ 
2008[ ). Data points are color-coded based on the galaxy’s 
^ value. The mean gas-to-dust ratio seen in our sample is 
^ 170, within the range found f or solar metallicity nearby 
galaxies ( Sandstrom et al. 2Q13[ ) and the Key Insights on 
Nea rby Galaxies: A Far-Inirar ed Survey with Herschel sam¬ 
ple ( Remy-Ruyer et al. [2014 ). The HCG galaxies show a 
relationship between their dust mass and gas mass that both 
matches the standard value and has little-to-no dependence 
on ^ or classification of MOHEG. This seems to indicate that 
the enhanced ^ values are due in part to a physical mecha¬ 
nism within the molecular gas, rather than an issue with con¬ 
verting CO luminosities into molecular gas masses. Addi¬ 
tional observations in ^^CO isotopologues (e.g. ^^CO and 
C^^O), as well as dense gas (e.g. HCN, HCO+, and CS), are 
necessary to confirm whether these systems require a different 
conversion factor, but the consistent gas-to-dust ratio appears 
to support a standard conversion factor. 


4. discussion 


4.1. AGNs in the HCG sample 

Table[^lists the properties of each of the CO-imaged galax¬ 
ies in our sample, including optical and radio signatures of 
AGNs. The detection of unresolved significant 3mm con¬ 
tinuum in 43% of the CARMA-imaged HCGs indicates the 
presence of AGNs, as the SFRs in these HCGs would not pro¬ 
duce sufficient free-free emission at 3mm continuum to be 
detectable by CARMA. Of the 7 HCGs detected in 3mm con¬ 
tinuum, 6 have 1.4 GHz radio detections as well (see Table [T] 
for fluxes and sources). HCG 82b is the only object detected 
in 3mm continuum but not 1.4 GHz. The optical nuclear spec- 
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Figure 4. The S^.gfxm/Ss.Gfxm versus Ss. 0 um/S 4 . 5 urn Spitzer IRAC col- 


2013 

1 com; 

et al. 

2009| 


comp ared to the Spitzer local voluW legacy sample (gray dots; |Dale| 
, with the so-called La cy wedge (circurtiscribing the A GN region 


.12004! [2007] T his plot 

onnson et al.| ( |200 /) and |Waik^ 


iws the g ai 

|et al.lpOlOj, highlighted as a yellow and black dashed liiie. Ihe CO-'imaged' 
HCCis are overplotted with the color of the point representing the level of 
SF suppression, the size of the points representing the mass of the molec¬ 
ular reservoir in each HC G galaxy, and the sh apes indicating whether the 
galaxy is a MOHEG from|Cluver et al.U2013J. Amongst the labeled HCG 
galaxies, squares represent MUHEGs (I/H 2 ,warm/^ 7 . 7 ^ 6 mPAH > 0.04) and 
circles non-MOHEGs. HCGs with the strongest suppression tend to be those 
that occupy the infrared gap in Spitzer color space. 

tra of the 3mm-detected HCGs all show signatures of AGNs 
(Table [^, either in the form of Seyfert-like spectra, low ion¬ 
ization nuclear emission line-region (LINER), or composite 


spectra ( [Martinez et al.||2010[ jCluver et al .J|20 13 j ) . However, 
in many of these cases, slow and fast shocks might be mim¬ 


icking these nuclear line ratios ( [Allen et al.|2008[ [Rich et ak 


[201 1[ Cales et al. 2015, submitted). 

There is an abundance of dynamically-excited features 
amongst our sample of warm H 2 bright HCG galaxies, es¬ 
pecially amongst galaxies with large values of ^ (although 
HCG 91a seems to also show signs of a broad blue line wing 
in its spectrum), which might be a sign that some of the AGNs 
seen are interacting with the nuclear molecular gas, possibly 
adding to the suppression (as is seen in NGC 1266; [Alatalo 


et al.[[20l3b| ). Deeper analyses of these systems might re 
veal molecular outflows, as evidence is mounting that out¬ 
flows are a commo n feature in objects with molecular reser¬ 
voirs near a n AGN ( [Garcia-Burillo et al.|2014|[Garcia-Burillo[ 
et al.[[2015[ ). A detailed analysis to determine whether these 


line wings and central broad lines are indeed due to molecular 
outflows would require higher resolution, higher sensitivity 
observations and is therefore beyond the scope of this paper. 


4.2. ^ and its connection to turbulence 

[Martig et al.[ ( [2Q13[ ) used simulations to predict that in a 
set of early-type galaxies with resolved molecular gas and 
SF, it was possible that massive bulges served to stabilize 
the molecular disks in these systems against gravitational col¬ 
lapse, reducing the SF efficiency of the molecular gas and 



8.5 9.0 9.5 10.0 10.5 11.0 11.5 12 

logCMgj [M®] 


Figure 5. The galaxy mass - star formation relation of a large sample of 
SDSS gala xies (gray contours), with SFRs and stella r masses derived using 
MAGPHYS jd^ Oanha et al.|2008||Chang et al.|2015l, showing the so-called 
“star formation mam sequence'' j Wuyts et al.] MTr Early-type galaxies are 
seen as the tail of high mass, low SEK. I he CAKMA-imaged HCGs are also 
overlaid, with colors coded based on their SF suppressions and symbol sizes 
based on their molecular gas masses. It is clear from this plot that the HCG 
galaxies with the most SF suppression are also the ones that are found to be 
farthest from the star formation main sequence. 


thus suppressing SF. For our sample, ^ does not correlate 
with the stellar mass of the system, the molecular gas frac¬ 
tion, (Mmol. gas/Mstar), or the visually-derived morpholog¬ 
ical type of the galaxy, suggesting that an alternative driver 
of suppression must be at play in these systems. HCGs that 
fall off of the K-S relation in Figure show similar behav¬ 
ior in Figure |^, suggesting that gravitational shears are not 
the main driver of the observed supp ression in our samp le (as 
was the case in early-type galaxies; Davis et ^[2014[ ). Al¬ 
ternatively, the derived properties for these objects might be 
inaccurate. We could be underestimating fpp, for instance, be¬ 
cause the molecular gas is not distributed in standard (Milky 
Way-like) molecular clouds, or because we have underesti¬ 
mated the molecular gas velocity dispersion (given the gas 
disruption that is common in HCGs). In fact, in these sys- 
te ms, it is possible that b oth diverge from the norm. 

Johnson et al.[ ( [2007 [ ) noted that a fraction of galaxies 


within HCGs seem to undergo a rapid transition between star¬ 
forming and quiescence, with a notable lack of galaxies with 
intermediate infrared colors. Figure [4| plots th e Spitzer In¬ 
frared Array Camera (IRAC; Fazio et al.|2004[) colors of th e 
CO-imaged HCGs originally plotted in [Lacy et al.[ ( [2004[ ), 
over laid with results from the Spitzer local volu me legacy sur¬ 
vey ([Dale et al.[200^ , as well as the HCGs from [Bitsakis et al. 


( [20111 ). Figure [4| shows that those of our galaxies located in 
the Spitzer IR gap tend to also have a high degree of SF sup- 
pression. 


Cluver et aL] ( [2013[ ) showed that the HCG galaxies that lie 
within the gap (and thus are rapidly transitioning; [Walker[ 


et al.pOlQ l also tend to be those with prominent warm H 2 


signatures that required mechanisms in addition to photon- 
dominated regions. Objects with ^ ^ 10 might be a ble to be 
reconciled with normal SF of Krumholz et al. ( [2012[ ) by revis- 
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ing our estimate of cTgas upward in the free-fall time equation, 
given that observatio ns of high-z galaxie s show Ogas with val¬ 
ues up to 50 km s“^ ICrescietaLpOO^. Tablepfshows what 


each galaxy requires to be reconciled with a normal SF ef¬ 
ficiency. Shocks heat the gas, injecting turbulence into the 
system, and boosting the molecular gas velo city dispersion 
by factors of a few to an order of magnitude (|Guillard et al.| 
[200^ . HCG25b, 40c, 57a, 79a, and 82b have ^ > 10, too 
high to reconcile with normal SF by applying a much higher 

<Tgas- 

SF requires gravitational binding energy to be greater than 
turbulent and radiative energies. If additional energy is in¬ 
troduced into the system, this balance is disrupte d, leading to 
SF becoming inefficient ( [Krumholz et al.|2012| ). One exam¬ 
ple of turbulence-induced SF suppression in the Milky way 
is the galactic center cloud GO.253-^0.016. This o bject shows 
evidence fo r a recent collision with another cloud dLongmor'e| 
et al.|2012| ) and has a lower SFE compared to similar objects, 
especially in regions with high velocity dispersion ( [Kauff-| 


mann et al.|2013 ). The shock in Stephan’s Quintet is dkected 
in CQ(l-O), but appears to have extremely weak assoc iated 


SF (|Konstantopoulos et al.|2014|). |Guillard et al.M2012a|) sug 
gest that ^ in this region could be a factor of 75 or higher, < 


result of the turbulence that has been generated by the shocks 
from the colliding galaxies. In fact, an increasing number 
of objects with strong turbulence have been shown to exhibit 


host NGC 1266 ([Alatalo et al. 

2015b[), and the radio galaxy 

3C326N (iGuillard et al.^Ols' 

). In this sample of objects. 


motions are observed associated with the suppression. The 
suppressed HCG galaxies in our sample seem to share some 
similarities with these systems, since many of them also con¬ 
tain gas exhibiting peculiar motions. 

How much energy injection is be required in order to sup¬ 
press star formation? If we assume that the molecular gas 
in the HCG galaxies are mainly rota tionally suppor ted, we 
can use the Toomre criterion (Q > 1; |Toomre|ri964| ), which 
describes the balance between rotation and turbulence, and 
gravitational binding energy, to determine the required en¬ 
ergy budget necessary to stabilize the molecular gas against 
collapse, effectively inhibiting star formation. To derive the 
necessary energy for s tability ((5=1), we use equation (7) in 
[Krumholz et al.| ( |2012| ): 


Q = 


y^2(y5 + 1) aQ 

TrfrSmol 


where a is the global gas velocity dispersion, Q is the rota¬ 
tion frequency, G is the gravitational constant, Smoi is the 
molecular gas surface density, and = 0 for the fiat part of the 
rotation curve. We then use force balance to determine the 
rotation frequency, ft. If we assume that all dis ks have equiv 
alent scale heights to the Milky Way (300 pc; |Burton|p"97T 


Malhotra|1995 ), the rotation frequency, ft, simplifies to: 




where H is the disk scale height. In order for the disk to 
be stable against collapse (Q= 1), assuming fiat rotation, the 
global molecular gas velocity dispersion must be: 


km / s 


7.0 


-'mol 


1 Mq pc 


-2 


1/2 


H 


300 pc 


1/2 


( 3 ) 


Assuming that the injected turbulent energy is 

^disk^ 

therefore be: 


^turb = the required Eturb must 


^turb ^ g ^ 20^^ /^^disk^ 


ergs 


H 


\3kpc/ \300pcy \lM 0 pc 


-'mol 


.-2 


( 4 ) 

And the turbulent luminosity (energy injection rate per rota¬ 
tional period) is: 


-^turb 

ergs/s 


: 4.8x10^^ 


f—V 

V300pcy 


f Rdisk \' 
\3kpc ) 


^mol 


IMo/pc^ 


( 5 ) 

Table lists these values for each of the HCG galaxies in 
our sam^e. Overall, the required gas velocity dispersions to 
stabilize the disk of the suppressed systems are reasonable, 
with the total injected turbulent energy representing ^1% of 
the total mole cular kinetic energy ( similar to what is seen in 
radio galaxies ; [Guillard et al.|2012b| ). Comparing the required 


turbulent injection to stabilize the disk to the extrapolated 


warm H 2 luminosity fromjCluver et al.|(|2013|), we find that 
in most cases, I/H 2 ,warm is within a factor of 3 of the required 
turbulent injection energy (save for HCG 96a and 96d, which 
is a known luminous infrared galaxy and star-bursting system, 
which likely has a gravitationally unstable disk). Many of the 
largest TH 2 ,warm/^turb values being found in galaxies with 
the largest suppressions, though there is also uncertainty in 
the extrapolation from the Spitzer IRS footprint to the galaxy. 


The fact that the H 2 luminosit; 

y, which is likely driven by tur- 

bulence ([Api 

pleton et al.|2006| 

1 Guillard et al.|2009[ 2012a|bl 

Cluver et a . 

j2013l) is comparable to the turbulent luminosity 


the star formation suppression that we are seeing in these sys 
terns, especially given the large contribution we expec t from 
[C II] cooling (up to a factor of 2 larger than ,warm; 


Guil- 


lard et al.]|2015| ), which will be discussed in detail in an up¬ 
coming paper (Appleton et al., in preparation). 

Other mechanisms than turbulence could be responsible for 
the suppression of star formation in those groups, like alter¬ 
native gas heating or tidal disruption processes. The cosmic 
ray heating of the ISM may be increased in interacting galax¬ 
ies ( jScalo & Elmegreen|20Q4| ), but even in systems showing 
extreme H 2 line emission and star formation suppression like 
Stephan’s Quintet and other groups or galaxy interactions, this 
heating mechanism would require an extremely high cosmic 


ray flux (jGuillard et al.||2009| |2012al [Peterson et al.||2012 
Cluver et al.[[2U13 ). Alternatively, the tidal field induced by 


the galaxy interaction can be responsible for some expansion 
of the gas on the external regions of the merger, and lower 
the average ISM pressure ( [Struck] 1 9'99} [Palous[2005 [ Renaud 
et al.|2009| , which could reduce the star formation efficiency 
locally. However this effect, which depends on the geomet¬ 
rical configuration of the tidal field and the relative position 
of the galaxy with respect to this field, is difficult to quan¬ 
tify without a proper numerical simulation of the interaction, 
and is generally thought to globally increase the compressive 
mode of turbulence (e.g. Renaud et al.|2014 ). 
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Schawinski et al.|(|2014) , overlaid with the 
Bitsakis et al. most suppressed 


Figure 6. Early-type (red contours) and late-type (blue contours) galaxies from the Galaxy Zoo JLintott et al.|2008||Schawinski et al.|2014|| Alatalo et al.|2014b) 
are compared to the HCGs. The color of the points represent th e level of SF suppres sion, the size of the points correspond to the size of the molecular reservoir, 
and the shapes indicate whether the galaxy is a MOHEG from|Cluver et al.U2013). Squares represent MOHEGs and circles non-MOHEGs. (Top le ft): The 
stellar mass vs. WISE [4.61-[121jum colors, overlaid with indigo lines dehnmg the IRTZ. (Top right): The u-r vs. r 4.61-ri21jum sequence identified by|Alatalo| 
with HCGs overplotted. (Bottom): The stellar mass vs. u-r colors of the Galaxy Zoo galaxies fromp ’ ’ ' ' ^ ’ ■ 

segment defining the optical green valley. The HCGs that are overplotted had both optical and IR colors derived from 
galaxies are the most likely to be in the WISE IRTZ and on or near the optical red sequence. 

Major mergers share many properties with HCGs, includ¬ 
ing the presence of violent interactions. In particular, the vio¬ 
lent interactions that thes e galaxies are exp eriencing result in 
shocks and turbulent gas ( [Cresci et al.J2009| ). Rather than sup¬ 
pressing SF, most of thes^bjects exhibi t super-efficient (but 
short timescale) SF ([Sanders et al.| 19911 [Sanders & Mirabelj 


[T^ [U et al.[[Mr2t [Lanz et al.[[2013[ ). The reason HCGs, 
which share many of these properties, have been found with 
suppressed SF rather than enhanced SF is possibly the result 
of the different interaction timescales, larger scale environ¬ 
ments, and gas properties. The CO(l-O) imaged HCG galax¬ 
ies in our sample are not observed during a major merger, and 
they are gas-poor relative to ULIRGs. 

The observed HCG galaxies have molecular gas masses and 
surface densities an order of magnitud e smaller than what is 
typically seen in inter acting galaxies ([Downes & Solomm 


1998t[Iono et al.|2009| . The molecular gas m many ULIRGS 


[& Scovillel 199^ [Wilson et al.|2Q08t[Ueda et al.|2Q14[ ). This 
likely means that although turb ulence can be a dis ruptive 
force, it is also transient (-1 Myr; [Guillardetal.[2009[ ). With¬ 
out a continuous input of energy, cooling lines are able to dis¬ 
sipate the turbulent energy rapidly, allowing the molecular gas 
in ULIRGS to restart forming stars quickly, due to the higher 
gas densities and shorter free-fall times. 

The timescale over which the violent interactio n is tak¬ 
ing place in major mergers is also short 10^ yr; [Hopkins 
et al.[[200^ Lanz et aL][2014[ ) compared to the total time in 
which gravitational forces withi n the compact group impac t 
the individual galaxies (^3 Gyr; [Plauchu-Frayn et ^[2012[ ). 
In major mergers, the gravitational enc ounters taking place 
involve only the two me rging galaxies ( Toomre & Toomre 
19721 [Privon et al.[[2013[ ). The dynamics within groups are 


IS more compact compared to HCG galaxies as well ( [Bryant[ 


much more complicated. HCG galaxies are interacting with 
the intragroup me dium, and their low densi ty gas is being ram 
pressure stripped ( [Rasmussen et al.[[2008[ ). Unlike the coa- 
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lescence timescale in major mergers, galaxies in HCGs are 
undergoing sustained gravitational encounters with the other 
group members, in which collisions do not necessarily result 
in coalescence. This extends the timescale of turbulent injec¬ 
tion, possibly resulting in the differences observed between 
major mergers and HCGs. _ 

An in-depth analysis by |Alatalo et al. ( 2Q14a ) of one of the 
most suppressed systems, HCG57a, suggested that an ongo¬ 
ing shock from a recent direct collision with HCG57d has 
continuously pumped turbulence into the system. H 2 and far- 
IR cooling lines have been unable to completely dissipate the 
extra energy, leading to the observed SF suppression. Once 
the shock has traversed the system (and thus is no longer 
pumping energy into the molecular gas), H 2 , [C ll], and other 
far-IR cooling lines should efficiently cool the gas, allowing 
SF to return to normal efficiency. It is possible that the group 
environment will extend this timescale as the group members 
continue to interact, shocking and expelling the interstellar 
media and ultimately transforming into early-type galaxies. 
This would be consi stent with the compac t group evolution 
picture suggested by Bitsakis et al.| ( |2010| ) that older groups 
contain a higher fraction of early-type galaxies. 

While it appears we have identified a sample of galaxies 
where turbulent energy has suppressed SF, the exact driver of 
this turbulence remains undetermined. Herschel data of far-IR 
cooling lines such as [C ll] and [O l] will be able to advance 
our understanding on the interplay between the molecular gas, 
SF suppression, and cooling mechanisms, potentially provid¬ 
ing us with a deeper understanding of how turbulence can im¬ 
pact the way in which galaxies form stars (Appleton et al., in 
preparation). 

4.3. SF suppression in warm H 2 -bright HCGs: a connection 
to transition ? 

Figure plots the stellar mass - SF main se quence for 
normal star-forming galaxies ( |Wuyts et al.||2Qll| ), c ompiled 
from the MAGPHYS -derived SFRs from SPSS DR7 ( |Abaza-| 
[jian et al.[|2009[[Qiang et al.|2015j ), with our CO-imaged ob- 
jects overplotted and color-coded by In this phase space, 
the HCGs with the most SF suppression fall the farthest from 
the main sequence, suggesting they are transitioning. Many 
of the other (non-suppressed) galaxies sit on the relation, in¬ 
cluding MOHEGs. In fact, the location of objects on the main 
sequence seems to mirror the situation for Spitzer infrared col¬ 
ors (Fig|^, strongly suggesting that the most suppressed ob¬ 
jects are also the ones that are currently transitioning between 
the star-forming spiral and quiescent early-type populations. 

Figure plots the CARMA-imaged HCG galaxie s on the 
optical ana IR color-magnitude diagrams presented in |Alatalo] 
|et al.l j lSOMbl ). The HCGs we have studied span a large range 
of optical and IR colors, although they are generally more 
massive than the underlying Galaxy Zo o distributions ( |Lin-| 
tott et al.||2008[ [Schawinski et al.|[2Q14b . One would expect 
that the total molecular mass or molecular gas fraction might 
be the determining factor of the position of the HCGs on the 
color-magnitude diagrams, but that is not the case. HCG 40c 
and 57a both contain large reservoirs of molecular gas, while 
appearing in the optical green valley or red sequence (com¬ 
pared to HCG 91a and 96a with equivalent molecular mass, 
which are quite blue). In fact. Figureshows that the mass of 
the molecular reservoir does not determine the color of the 
galaxy. Spearman rank tests were used to search for cor¬ 
relations between ^ and other galaxy properties.The only 
galaxy property significantly correlated with ^ was color. 


Including HCG 68a, a Spearman correlation of -0.71 with a 
value = 0.0044 was found when investigating ^ and u-r 
color. If HCG 68a is removed, we find a y-(u-r) Spear¬ 
man correlation of -0.71 (equal to the test with HCG 68a 
included) with a /7-value = 0.0011. Without HCG 68a, 
([4.6]-[12]/im) also shows a Spearman correlation of 0.80 
with ap-value = 0.0067. With or without HCG 68a, the optical 
and IR colors were the only galaxy properties shown to signif¬ 
icantly correlate with y. The extreme IR colors of HCG 68a 
might be due to a buried AGN (Alatalo et al. 2015c, in prep). 

These results seem to suggest that a galaxy does not need 
to shed its IS M to then quench SF (the standard galaxy tran¬ 
sition picture; [Hopkins et al.|2006 ) to morphologically trans¬ 
form. Instead, changing the state of the molecular gas can 
act to quench SF before the ISM has been completely shed. 
It is thus possible that in many of our CO-imaged galaxies, 
the molecular gas has been rendered infertile due to shocks 
pumping turbulence into the system. The quenching of SF, 
and the beginning of the transition across the green valley, 
occurs before the galaxy loses the majorit y of its molecular 
ISM. This is consistent with the findings of [Leon et al.|(|1998|) 
and Martinez-Badenes et al.j ( |20I2| ) that HCG galaxies con- 
tain comparable molecular gas reservoirs to is o lated galaxies 

! Lisenfeld et alJ20lT| Saintonge et al.|2QlT^ . [Alat^o et al.j 
2014b| ) posited that the IRTZ was a manifestation of evo- 
ution, representing the stage in which a galaxy is actively 
shedding its ISM, and that this phase follows SF quenching 
(traced by the u-r colors). In this scenario, first the galaxies 
move toward the elbow in Fig [^, where SF is suppressed, 
and then they move into the IRTZ, where they quench S F and 
move to the optical red sequence. The Hi findings of jSerra j 
et al.| ( |20f^ of that 40% of field early-type galaxies still con- 


tam non-negligible reservoirs of neutral gas. New work on 
poststarburst galaxies found that many poststarburst galax¬ 
ies still contain non-negligible m olecular reservoirs ( [French 
et al. [2015} [Rowlands et al.|2015[ ), confirming the suggestion 
that SF quenching can take place before galaxies expel their 
molecular interstellar media. It is currently unclear how of¬ 
ten galaxies transition in this fashion, but further observations 
of IRTZ and poststarburst galaxies, combined with our new 
observations of at least some HCG galaxies (certainly those 
that are warm H 2 -bright and CO-bright), have shown that this 
“quenching first” path contributes to the population of transi¬ 
tioning galaxies. 

The presence of suppressed objects within the IRTZ points 
to a possible new method for identifying galaxies most likely 
to be in a phase of inefficient SF, before the exhaustion of 
their molecular gas. By directly imaging the CO in molec¬ 
ular gas-rich galaxies that also appear in the IRTZ, we can 
pinpoint a population of galaxies likely to exhibit SF suppres¬ 
sion. This result suggests that future studies might identify 
larger samples of suppressed galaxies by selecting galaxies 
based on their optical and IR colors. A larger sample of sup¬ 
pressed systems thus allows us to study how the injection of 
turbulence n ot only impacts the e nergy balance that dictates 
our SF laws ( [Krumholz et al.[2QT^ , but also allows us to study 
how the neutral ISM is exhausted as a galaxy transforms from 
a spiral into an early-type galaxy. 

5. SUMMARY 

• We have used CARMA to map CO(l-O) of galaxies 
within 12 HCGs, many with elevated warm H 2 emis¬ 
sion, detecting molecular gas in 14 galaxies, and unre¬ 
solved 3mm radio continuum in 7 (consistent with the 
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presence of AGNs). Figures [CT- C12 show the molec¬ 
ular gas data for each galaxy, including the momentO 
map, moment 1 map, channel maps, PVDs, and inte- 
grate d spectra. A compar ison to the single-dish data 
from |Lisenfeld et alT] ( |2014| ) shows that our observations 
have not resolved out large fractions of the flux. 

The molecular gas morphologies of our HCG ga laxies 


Alat- 


(using the metr ic set for the ATLAS galaxies in^ _ 

alo et al.|2Q13| ) indicate that HCGs are consistent with 
the distribution of gas morphologies foun d in early-type 
and spiral galaxies rather than ULIRGs ([Wilson et al. 
[2008i|Uedretal.|2014| i. 

• We have shown that a large proportion of our CO- 
imaged HCG galaxies exhibit SF suppression (^) 
when plotted relative to both the K-S relation (|Ken-| 


[ nicutt[ 1 

|I998[) as 

well 

Krumholz et al.[ 

(20T2| 

this sample is (J 



The mean SF suppression for 
I0±5, and exhibits a bimodality. 
The most extreme objects (HCG 25b, 40c, 57a, 79a & 
82b) exhibit ^ > 10, and have molecular gas depletion 
timescales t^ep ^ lOGyr. 


• The mean gas-to-dust ratio for the CO-imaged HCGs 
is a round 170, within the range found in normal galax - 
ies ( [Sandstrom et al.|[2013[ |Remy-Ruyer et al.||2014 ). 
We do not believe that the observed ^ is due to an 
incorrect Lco-M(Y{ 2 ) conversion factor (which would 
appear as a highly discrepant gas-to-dust ratio in sup¬ 
pressed galaxies). 


• A non-negligible fraction of our CO-imaged HCG 
galaxies contain ed substantial warm H 2 emission ( |Clu 
ver et al.|2013| ), consistent with there being shocks m 


jecting substantial turbulence into these systems, and 
the turbulent energy required to stabilize the molecu¬ 
lar gas against collapse appears to agree within an or¬ 
der of magnitude with the warm H 2 luminosity. As 
has been seen in the Milky Way ([Kauffmann et al. 
[ 20 T 3 ]), NGC 1266 dAlatalo et al.|2015bD , and 3C326N 
( [Guillard et al.||20i5| ), the additional turbulence could 
upset the energy balan ce that dictates the rate of SF 
( Krumholz et al.|20T2 ), thereby suppressing SF. 
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bars arise during gravitational enc ounters ( [Athanassoula 


1996t lAthanassoula & Bureau||19^ , it is unsurprising that 
■Oc is morphologically classed as a bar+ring. 


A. DETERMINING THE AGN CONTRIBUTIONS WITHIN 
CO(l-O) IMAGED HCG GALAXIES 

In Table we presented the AGN classifications of the 
galaxies in our sample, based both in optical emission line 
diagnostics as well as radio continuum. Such sources can 
significantly contribute to the infrared emission of the galax¬ 
ies, leading to an overestimation of their star formation rates, 
if they are not accounted for. To disentangle the fraction of 
AGN contribution in the total infrared luminosity we have fit¬ 
ted the galaxy observed infrared (8-500/im; purely dust emis¬ 
sion) spectral en ergy distributions with DecompIR dMul- 
[laney et ar]|2011| ). This code simply fits the observed fluxes 
with sets of host-galaxy AGN component templates, and 
estimates the contribution of the AGN to the total infrared lu¬ 
minosity. The AGN templates are described by broken power- 
laws at around 40/im that fall steeply above that. From this 
analysis we find that only three of our sources have signifi¬ 
cant AGN contribution at the infrared bands, HCG 68a with 
9%, and HCG 91a with 35%. 


B. COMMENTS ON INDIVIDUAL GALAXIES 


HCG25b: Figure [Cl shows that HCG25b is an edge-on 
galaxy. Deep optical imaging by Eigenthaler et al.j ( |2015| ) 
also showed that HCG 25b is interacting strongly with 
HCG 25f. A tidal tail connects the two, and the polyaromatic 
hydrocarbon (PAH) emission shows a small tida l feature. 
HCG25b is also a MOHEG ( jCluver et al.||2013] ) and is a 
transitioning galaxy in the IR (i.e., it lies in the optical red 
sequence as well as the Spitzer IR gap and the WISE IRTZ). 
HCG 25b also likely contains a n AGN, classified thr ough 
optical emission line diagnostics ( |Martmez_et al.||2Q10| ), and 
the presence of 1.4 GHz nuclear emissioij^ The molecular 
gas is morphologically classed as a disk. 


HCG 40c: Figure [C2] shows that HCG 40c is an edge-on 
galaxy strongly detected with CARMA. While HCG 40c lies 
in the red part of the WISE IRTZ, it is found in the optical 
green valley, though near the red sequ ence. It also has a radio 
core and optically identified AGN (Martinez et al. |2010| ). 
Two Micron All-Sky Survey (2MASS) imaging a lso shows 


et al 


that HCG 40c is tidally interacting with HCG40e (Skrutskie 
:][2006| ). The PVD of HCG 40c (Fig|g also shows a 
ficant bar (seen as the large velocity structure with very 


sigmticant bar (seen as the large velocity structure w itb very 
little positio n shift), similar to what was seen by [Alatalo 


et al. 


( 2013| ) in several ATLASearly-type galaxies. This is 


likely what is responsible for the slight appearance of broad 
wings in the CO spectrum. Given the frequency with which 


^ Nuclear 1.4 GHz emission can also be due to star formation jCondonj 
Thus, detecting 1.4 GHz emission in these objects doe s not con- 


hrm t he presence of an AGN without morphological confirmation JBest et al.| 
|2005| and thus are just suggestive of their presence. 


HCG 47a: Figure shows that HCG 47a is an oblong 
spiral (in both PAH and optical images from SDSS; [York 


et al. 


_ that is tidally interacting with HCG 47b. 
HCG 47a is also on the cusp of being considered a M OHEG 
(rH2,warm/i7.7^m PAH = 0.035; [Cluver et aT]|2013| >. The 
center, while vo id of molecular gas, contains a spectrally 
classified AGN ( [Stern & Laor||201^ as well as 1.4GHz 
emission. The CO(l-O) emission traces the oblong spiral, 
and so the molecular gas is morphologically classified as both 
a ring and a spiral. 


HCG 55c (Fig |C4[ ) is part of the chain of galaxies 
VV 172/Arp 23S[^ HCG 55 is the most distant HCG 
that we have imaged with CARMA. The emission in this 
galaxy is resolved and spans approximately 3 beam widths. 
HCG 55c was not detected in 1.4 GHz or 3mm continuum, 
and does not have a published optical spectrum (to determine 
the presence of an AGN). However, HCG 55c is at the edge 
of the Spitzer gap, and is in both the optical green valley 
and the WISE IRTZ. The molecular gas in HCG 55c is 
morphologically classified as a disk. 


HCG 57a is found to have some of the most complex molec¬ 
ular gas dynamics in the sample, including three distinct 
kinematic components (Fig reffig:hcg57), and is interacting 
with HCG 57d. HCG 57a is also a MOHEG dCluver et al.] 
20f^ , and an in-depth dis cussion of both HCG 5 7a and 
HCG57d can be found in [Alatalo et al.[ ( [2014a[ ). The 
molecular gas in HCG 57a is morphologically classified as 
mildly disrupted. 

HCG 57d (Fig [C5[ ) was not within the Spitzer IRS footprint, 
and thus we do not have warm H 2 information on this source. 
An in-depth discussion of both HC G 57a and HCG57d can 
be found in Alatalo et al.l ( 2014a[ ). The molecular gas in 
HCG 57d is morphologically classified as a ring. 


HCG 68a (Fig [C6[ ) was previous ly un-detected by the 
IRAM 30m in the ATLAS^^ survey (Young et al. 201 1[ ), but 
was later detected by [Lisenfeld et al.[ ( [2014[ ). The CARMA 
observation helps shed light on why this was the case. The 
molecular gas in HCG 68a is not only compact (unresolved 
by the CARMA beam), but also very broad (Ai;^2000 
km s“^). HCG 68a also appears to have rotation velocities in 
its CO e mission that makes it an outlier on the Mgai-Prot,co 
relation ( [Davis et ~aL [2011[ ). HCG 68a and HCG 68b were 
also both detected in 3 mm continuum emission. DecompIR 
( [Mullaney et al.[[MTT[ ) suggests that about 9% of the far-IR 
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It contains a galaxy with a highly discrepant redshift (HCG55e, 


r; = 36880 km s 


Sargent 


19681. 
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emission (|Bitsakis et aL||20141) originates from the AGN in 
this system. The AGN in HCG 68a has an 2-10 keV X-ray 
luminosity of Lx = 1.6 x 10^° ergs s“^ ( Evans et al.|2010| , 
which could account for the majority of the far-IR emission 
if the obscuring column of molecular gas is sufficiently 
high (which the CARMA observations suggest might be the 
case). HCG 68a is also one of the strongest MOHEGs, with 
^H2,warm/^7.7Aim PAH = 0-741, and is an outlier in IR color 
space amongst the CO-imaged HCG galaxies. HCG 68a 
has Spitzer colors corresponding to an early-type galaxy, 
though these colors might also be due to a buried AGN with 
very few intermediate aged stars (Alatalo et al. 2015c, in 
prep). HCG 68a is also the only galaxy that is found with 
WISE colors complet ely consistent with the elliptical sample 
( [Alatalo et al.|2014b| ), but sits within the optical green valley. 
The molecular gas in HCG 68a is morphologically classified 
as a disk. 


HCG 79a (Eig |C7|), a member of Seyfert’s Sextet ( jSeyfertj 


1951p is a near edge-on red early-type galaxy with a promi 


nent dust lane. HCG 79a is optically classified as having an 
AGN and is found in the WISE IRTZ. The molecular gas in 
HCG 79a is morphologically classified as a disk. 


HCG 82b (Eig |C8| ) is the only CO-imaged HCG that was 
successfully detected in 3mm continuum but not in 1.4 GHz 
continuum. Despite th is, there is an optically identified AGN 
( [Martinez et al.|20T0| , so the 3mm conti nuum is most likely 
due to an AGN HCG 82b is a MOHEG ( [Cluver et al.[[2013| ) 
located in the red seq uenc e, as well as the Spitzer gap and 
the WISE IRTZ. Eigure JCS] shows HCG 82b contains a stellar 
bar, where the molecular gas does not appear aligned with 
the bar. The lip seen at the edge of the molecular disk seems 
consistent with a warp, and thus the molecular gas in this 
galaxy is morphologically classified as mildly disrupted. 


HCG 91a (Eig [C9[ ) is a nearly face-on spiral galaxy that 

is blue in optical and IR colors, with the CO(l-O) being 

brightest in the northeast quadran t of the galaxy. Deep Ha 

imaging from [Eigenthaler et al.[ ( [2015| ) seems to indicate 

that this might be the location of a bow shock. I t is the 

only galaxy in this survey classified as a Seyfert 1 ([Cluver 
^ .... 


et al. [2013| ). The deep optical imaging of HCG9ia also 


indicates t hat it is undergoing a sig nificant interaction with 
HCG 91c ( [Eigenthaler et al.[ 2015j ), showing multiple tidal 
tails, including on e connecting the two interacting galaxies. 
Vogt et ah] ( [2015] ) also showed that the CO and optical line 
velocities were offset from that of the Hi. HCG91a has 
many different kinematic components but appears to have 
a blue-shifted line wing that runs directly north-to-south 
from the nucleus of the galaxy. DecompIR ( [Mullaney et al. 


2011 ) suggests that up to 35% of the far-IR pitsakis et ^ 


2014) from the center is due to an AGN. The molecular gas 


m HCG 91a is morphologically classified as a spiral. 


HCG 95c (Eig [C10[ ) mainly shows regular rotation, but 
the CO(l-O) emission extends into the tidal tail between 
HCG 95a and 95c, which is also seen prominentl y in t he 8/im 
non-stellar emission in the channel map of Eig [C10[ Shells 
are also present in the optical light. Spitzer IRAC colors 
place this galaxy on the edge of the Spitzer IR gap closer to 
those of dusty spirals, but its H 2 /PAH ratio indicates that it 

although HCG 79c later identified to be a background galaxy 


is a MOHEG ( [Cluver et al.[[2013j ). Unlike other transition 
galaxies, this one has a high specific star formation rate. 
HCG 95c co ntains radio emission, a nd is spectrally classified 
as an AGN ([Martinez et al.[[2010[). The molecular gas in 
HCG 95c is morphologically classified as mildly disrupted. 


HCG 96a (Eig [Cll[ ) is the only CO-imaged HCG galaxy 
that could be iden tified as containing a n AGN solely from its 
mid-IR spectrum ( [Cluver et al.[[2013] ) , that is also classified 
as an AGN by optical spectroscopy ([Martinez et al. 2010[ ). 
Deep optical imaging ( [Eigenthaler et al. [2Q15[ ) shows that 
HCG 96a is undergoing an interaction with HCG 96c, with 
multiple tidal tail s present, including stellar light c onnecting 
HCG 96a & 96c ( [Verdes-Montenegro et al.[[T997] ). Optical 
and WISE IR colors are all consistent with a star-forming 
galaxy, thoug h this object is lo cated in the A GN wedge in 
Spitzer colors ( Lacy et ar|[2004[ ). Eigure [CTT[ shows that the 
molecular gas contains multiple components, including spiral 
structure, a bar component and a ring. Thus the molecular 
gas in HCG 96a is morphologically classified as both a spiral 
and a bar-^ring. 


HCG96c (Eig [Cll[ ) is interacting with HCG96a, as men¬ 
tioned above, and appears in the optical green valley, and 
has Spitzer and WISE colors consistent with star-forming 
galaxies. The 1.4 GHz emission indicates that HCG 96c 
might contain a radio-bright AGN. HCG 96c also has a broad 
CO(l-O) spectrum, but that is possibly due to its edge-on ori¬ 
entation. The molecular gas in HCG 96c is morphologically 
classed as a disk. 


HCG 100a (Eig C\2) exhibits prominent spiral structure, 
with prolific SE activity in its center (seen in optical light 
and PAH emission in Eig[C^. HCG 100a does not exhibit 
any outward signs of interaction, such as tidal tails, and 
is identified as star-forming in both optical and IR colors. 
HCG 100a is spectra lly identified as containing an AGN 
( [Martinez et al.[[MT0| ), consistent with its detection in 3mm 
continuum and 1.4 GHz emission. The majority of the 
molecular gas in HCG 100a seems to be undergoing regular 
rotation, but there also seems to be a minor-axis component 
that includes higher velocities (seen as wings in the CO(l-O) 
spectrum), possibly due to a bar or an outfiow. Despite the 
presence of a putative minor axis component of the molecular 
gas in HCG 100a, we morphologically classify this molecular 
gas as a disk. 


C. COMPLETE FIGURES FOR EACH HCG 

Eigures for each individual HCG are shown in Eigures 
Cl - [C12 Eor each HCG, we show five figures, as well as 


an additional figures for HCGs where CARMA detected 2 
galaxies: HCG57d (Eig[CT^) and 96c (Eigs. [CT3| ) &[CT^. 
The outlined box in the image (white dotted line) represents 
the area of the corresponding momentl map. North is up and 
East is left in all images. 


Top left: The CO(l-O) integrated intensity (momentO) map 
(white contours) overlaid on an optical 3-color image, either 
g,r,i from the Sloan Digital Sky Survey (HCG 25, HCG 47, 
HCG 57, HCG 68, HCG 79, HCG 82, HCG 95, HCG 96, 
HCG 100), the Digitized Sky Survey (HCG 40, HCG 55), or 
the Swift archive (HCG 91). 
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Top right: The CO(l-O) mean velocity (moment 1) map, 
overlaid with the momentO map (white contours). The 
velocities are relative to the systemic velocity of each source. 

Middle left: The integrated CO(l-O) spectrum, which was 
constructed by summing all pixels in each channel using the 
momentO mask as a clip mask. 

Middle center & right: The position-velocity diagram 
(PVD), taken by slicing the CO(l-O) cube in a plane and 
summing over a specific region. The momentO map (middle 
center) outlines the PVD integration area (dashed black line 


for the center, dotted gray lines for the boundaries), and the 
corresponding PVD shows the velocity structure tangent to 
the velocity slice (middle right). 

Bottom: The CO(l-O) channel maps are overlaid on the 
Spitzer 8/im nonstellar emission. The Spitzer 8/im nonstellar 
maps were created by subtracting a scaled 3.6/im map from 
IRAC from the 8.0/im IRA C map using the s cale factor of 
0.232 for late-type galaxies ( |Helou et al.j2004| ). The shading 
corresponds to a signal-to-noise ratio of 3 in each channel 
(except where noted otherwise), with additional contours 
in either Icr or 3cr steps. The panel that corresponds to the 
systemic velocity has its velocity labeled in red. 
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Figure Cl. HCG25. Channel map contours are in Icr steps. Elements of this figure are described in 
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Figure C2. HCG40. Channel map contours are in Zg steps. 
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Figure C3. HCG 47. Channel map contours are in Scr steps. 
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Figure C4. HCG55. Channel map contours are in la steps. 
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Figure C5. HCG 57. Channel map contours start at ±2.5(7 and are in 1(7 steps. The PVD of 57d can be found in Fig |C13| 
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Figure C6. HCG68. Channel map contours are in Icr steps. 
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Figure C7. HCG79. Channel map contours are in 3cr steps. 
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Figure C8. HCG 82. Channel map contours are in 3cr steps. 
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Figure C9. HCG91. Channel map contours are in Str steps. 
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Figure CIO. HCG95. Channel map contours are in 3a steps. 
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Figure Cll. HCG96. Channel map contours are in 3a steps. Channel maps and a PVD of 96c can be found in Figures [cT4] and [cT3] respectively. 
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Figure C12. HCG 100. Channel map contours are in Zg steps. 
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Figure C13. PVDs of HC G 57d (top) and HCG 96c (bottom). The molecular gas in HCG57d appears to be consistent with a ring (seen in PAH emission in 
Figure [C^ and described in |Alatalo et al.|201^ . The molecular gas in HCG 96c seems to be consistent with a rotating disk. 
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Figure C14. HCG 96c. Channel map contours are in la steps. 







































